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I. INTRODUCTION 



Flavor Changing Neutral Currents (FCNC) in general, and of the top quark in particular, 
play an important role as a testing ground for the Standard Model (SM) and for New Physics 
(NP). Because of its large mass, the top quark can decay into all other quarks, accompanied 
by gauge or Higgs bosons, as well as into new particles predicted by NP models. The interest 
in FCNC for top quark physics stems from the facts that: 

1. The scale of NP is closer to the top quark mass more than to any other quark. 



2. In the SM, the FCNC two-body processes t — > c</,7, Z, i7 are absent at tree- level 
and are highly suppressed by the GIM mechanism at one loop. Their branching 
ratios predicted in the SM are of the order of 10~ n to 10~ 14 [l], Q]> f ar away from 
present and even future reaches of either the Large Hadron Collider (LHC) 0, 0] or 
the International Linear Collider (ILC) There are many models of NP in which 
the branching ratios for the above 2-body FCNC decays are much larger than those 
obtained in the SM (see e.g. jflH] and references therein). 

In addition to the two-body rare decays of top quark, some of its rare three-body decays 



e.g., t —>■ c WW x cZZ, bWZ have been considered in the literature within the SM 



and for NP 



id flll | . 



decays in the SM but some of them get comparably large within models of NP, such as 
two-Higgs-Doublet jllj . especially after including finite-width effects jloT ]. 

In this paper we will analyze another three-body rare decay, namely t — > egg within the 
SM framework and compare it to both t —>■ eg and t — > cqq, q = u. The main motivation 
for such a calculation comes from the phenomenon in which higher order dominates over a 
lower order rate, as observed in the c, b and in other systems. 

For the case of charm decays, the short distance contribution to c — > wy, exhibits a huge 
enhancement over the lowest order penguin diagrams 12j. One can argue then that even 
higher order short distance are not that important. Unfortunately for radiative D decays, 
even this enhancement is overshadowed by much larger long distance terms. 

For 6-quark decays, a study of the next-to-leading logarithmic (NLL) QCD corrections 
[l3l |. yielded Br NLL (6 — > sg) ~ 5.0 x 10~ 3 , whereas the leading-logarithmic (LL) result 
Br LL -> sg) w 2.2 x 1(T 3 , in spite of an a s /ir suppression. This large correction is 
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dynamical in nature, since it is due to a arge ratio of Wilson coefficients evaluated at m&, 
C 2 /C$ ~ 7 (for earlier references see jl4, 13 13)' Higher order dominance may also 
become substantial in the decays of a sequential fourth generation quark b' , if it exists 



15 



18 
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More recently |2Cj, the one-loop, three-body, rare top quark decay t — > u\U 2 u 2 , where 
Ui = u,c was calculated in the SM and found to dominate over the one-loop, two-body, rare 
t — > u\g decay, by about one order of magnitude, although the latter is of lower order in a s . 
Later on, we will comment about their result. 

The purpose of this article is to evaluate and discuss the higher order dominance issue in 
rare top decays among t — > cg,t — > cqq, and t — > egg. The present calculations are within 
the SM and while, as discussed above, t — > eg and t — > cqq, q = u were calculated before, 
this is the first calculation of t — > egg in the SM. 

The remainder of the paper is organized as follows: In Section II we present the calculation 
of t — > egg, in Section III the decay t — > cqq, q = u is evaluated, and in Section IV we 
conclude. The Appendix includes the one-loop functions which appear in Section II. 



II. CALCULATION OF t -> egg 

The one- loop t — c — g*(k) and, in general, the oi — q 2 — g*{k) vertex function can be 
expressed, using Lorentz and gauge invariance, as [2jJ 

T M = F 1 (k 2 ) (k 2 ^ - kj) P L - %F 2 {k 2 )m t <v k v P R , (2.1) 

where Pl,r = (1 T 7s)/2 and m c = are assumed. The functions Fi and F 2 are called 
charge-radius (or monopolar) and dipole moment (or dipolar) form factor, respectively. Note 
that this is not the most general vertex function. There are two more form factors, namely 
Fin (right-handed monopolar) and F 2 l (left-handed dipolar), which are both proportional 
to m c /m t so that for the sake of simplicity we omit them here (see [3 f° r the details). In 
our numerical analyses, however, all contributions are retained. 

While Fi(k 2 ) contribution to q\ — > q 2 g, q\ = t, q 2 = c (i.e. for a real gluon) vanishes, 
both Fi{k 2 ) and F 2 {k 2 ) give non-zero contribution to qi — > q 2 qq and to qi — > q 2 gg. Of 
course, the vertex functions for the three-body processes are more complicated than 
just above, but if F 2 < F±, there is a chance that the three-body modes will be of the 
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same order, or even larger than t — > eg. 



c c 




FIG. 1: The one-loop contributions to t — > egg in the SM in the 't Hooft - Feynman gauge. The 
ghost contributions are depicted in Fig. [21 The first 1-10 diagrams are the vertex diagrams, the 
diagrams 11-14 are the one-particle irreducible (1PI or box) diagrams, and the rest 15-30 are the 
t — c self energy diagrams. The crossed diagrams are also shown explicitly. 

In the SM, the decay t —* egg occurs at one-loop level and the Feynman diagrams 
contributing to the decay are depicted in Fig. \I\m the t' Hooft - Feynman gauge (£ = l). 1 
The G field in the diagrams is the unphysical part of the Higgs field. The polarization sum 

1 Note that there will be no cross term for diagrams with triple gluon vertex since it is already counted in 
the vertex factor. 
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of the gluons is naively 



5^e*0, X)e v (k, A) = -g^. 



(2.2) 



However, it is well known |22j, |23[ that in cases where there are two or more external gluons 
(either for tree, like gq — > gq, or for loop diagrams as in the present case), the above 
sum leads to violation of gauge invariance. The problem is alleviated either by choosing a 
transverse polarization sum, or by introducing ghost fields to get rid of the unphysical gluon 
polarizations while keeping the simple polarization sum above. The ghost contributions are 
shown in Fig. |21 




FIG. 2: The one-loop ghost contributions to t — > egg in the SM in the 't Hooft - Feynman gauge. 
There is a ghost contribution for each diagram of Fig. ^with triple gluon vertex. 



Fig. |21 requires some explanation. The ghosts couple only to a gluon, i.e. there is a 
(^-ghost-ghost coupling but there is no coupling to quarks. Therefore, there is a ghost 
anti-ghost diagram for each diagram with a triple gluon coupling. Furthermore, since the 
ghost is not its anti-particle there are six diagrams with ghosts as in Fig. 121 and there is no 
statistical factor after phase-space integration, unlike the case of two final gluons where a 
statistical factor of 1/2 is inserted following phase space integration. Note also that there 
are no interference terms between ghost and ghostless amplitudes. 2 

In our calculation we make the GIM mechanism manifest by dropping terms independent 
of internal quark mass. 3 This will substantially simplify the calculation. The divergent parts 
originate only from the unphysical Higgs loops. As the cancellation of ultraviolet divergences 
should happen at the amplitude squared level for t — > egg in the 't Hooft - Feynman gauge, 
the analytical output is so long that the ultraviolet finiteness has to be checked numerically. 



We carried out the calculation in D-dimension using dimensional regularization 



2J and 



2 The ghost holds, although they are bosons, obey Fermi-Dirac statistics and will have an overall -1 factor 
for each closed loop, similar to a fcrmion loop. 

3 This simplification applies to the divergent parts as well since the SM without the unitarity of the CKM 
matrix is non-rcnormalizable. 
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keep the decay width formulas as functions of e = 4 — D. Then we check the stability of 
the decay width by varying e. Of course, this method is limited by the achievable numerical 
precision of the program used. We have done our calculation with the softwares FeynArts and 



FormCalc 



261 ] and done some partial cross-checks with FeynCalc 27j]. With FormCalc 



the numerical check of ultraviolet finiteness is controlled by a parameter and thus easier to 
check. We have tested and confirmed the finiteness of our results for t — * egg (as well as for 
t — > cqq, q = u and t — »■ eg). 

In addition to ultraviolet divergences, we must deal with infrared and collinear divergences 
which exist in the calculation of t — > egg. There are three possible cases producing such 
singularities: a) the configuration of having one of the gluons travel parallel to the charm 
quark, b) the configuration of having the gluons traveling parallel to each other, or c) the 
fact that one of the gluons in the final state could be soft. While cases a) and b) are related 
collinear divergence, case c) is related to the infrared singularity. To cure case a) we take 
a non-zero mass for the charm quark in our numerical calculation even though we present 
some of our analytical results in the m c = limit, just for simplicity. 

There are two ways to approach cases b) and c). One can either do an extensive study by 
including QCD corrections consistently of the same order in the perturbation theory, this 
case up to the order of a 2 s requires including the interference terms from two-loop t —* eg 
(one-loop QCD corrections) with one-loop t — > egg; or one can exclude the part of the phase 
space containing singular points by simply imposing cuts on the kinematics of the process. 
The former method can be achieved by dimensionally regularizing the phase space integrals 
with D = 4 — eiR and making the divergence manifest. Then the cancellation of infrared 



H 29, Q- 



singularities is guaranteed by the Kinoshita-Lee-Nauenberg (KLN) theorem 
Basically one has to carry out a study for top quark decays similar to the b —>■ sg evaluation 
done by Greub and Liniger[l^]. We chose to proceed by cutting the "dangerous" integration 
limits as discussed below. Our calculation is therefore more in the spirit of 

When we carry out integration over phase space by using the momentum delta functions 
and azimuthal symmetry, we are left with only two non-trivial integration over variables, 
say, the energy of one of the gluons (E 3 = k®) and the energy of the c-quark (E c = k®), 
where the energies are in rest frame of the decaying top quark. Then in (E c , E 3 ) space, the 
region of integration becomes a triangular shape. The singularities discussed above lie at 
the boundaries of the region, infrared singularities at the vertices and collinear singularities 
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along the boundary lines. Thus we have to impose cuts on both E c and E 3 . We will discuss 
this point further later in the section. Note that since we constrain our phase space, our 
result for Br(£ — > egg) should be seen as approximate and more conservative. 

Let us now present some analytical expressions. For the sake of simplicity of the pre- 
sentation, the masses m c , and m s are assumed to be zero, though we included their 
contributions in our numerical study. From Figs. ^ and EJ the amplitude for the decay can 
be written in a compact form as 

aa s V tb V* b 



A 



self 



[-F 2 m t t l2 R 2 R 4 - 2F 1 (SP 4 - SP 5 )tR 6 R 4 



A 



4m^m^s 2 3sin 9 w tti 2 R\R 4 
+s 23 SP 1 R 1 (2F 9 R 5 - F 12 m t R 3 R 4 ) + F 3 t 12 R 7 + 2F 4 R, 
aa s V tb V* b 



vert 



.4 



ipi 



2m^s 23 sm 2 9 w tt 12 RiR 4 
—F 3 ti 2 Ri 4 + R 4 (m t s 23 ti 2 (2Fi 5 R w 
aa s V tb V* b 



2mL sin 2 6> 



^ghost 

where 



[2Fgi?i 7 

w 0111 u w 
— 2F\ 2 R 22 ) + 2F 4 R 23 

aa s V tb V; b 



[2^(^3^*12^(9) - F 9 s 23 SP 1 R n ) - 2F 4 R l3 
F 2 R l2 ) + 2F 1 tR 15 ) + F l2 m t R w } , 
F 3 R l9 + m t (-2F 15 R 18 + 2^3^ + F 2 R 21 



^m 2 t m\y sin 9 w t\ 2 



2F\R 24 ] , 
-R25 [-P9-R27 ~ 2F 12 m t R 2& ] , 



(2.3) 



Fx 


= u{k 2 


F 3 


= u(k 2 


Fr, 


= u(k 2 


F- 


= u(k 2 


F, 


= u(k 2 


F n 


= u(k 2 


F± 3 


= u(k 2 


F X5 


= u(k 2 



0)P R f(k 3 )u{k 1 ,m t ), F 2 = u{k 2 ,0)P R f{k 3 )f(k 4 )u{k 1 ,m t ), 
0)P R f(k 3 )f{h)^ 3 u(h, mt), F 4 = u(k 2 , 0)P R f(k 4 )u{k 1 , m t ), 
0)P L f(k 3 )u(k 1 ,m t ), F 6 = u(k 2 ,0)P L f(k 3 )f(k 4 )u(k 1 ,m t ), 
Q)P L f(k 3 )f(k 4 )ft 3 u(k 1 , m t ), F 8 = u{k 2 , 0)P L f(k 4 )u(k 1 ,m t ), 
0)P R l / k 3 u(k 1 , m t ), F w = u(k 2 , 0)P L u(k 1 , m t ), 
0)P L ] ! k 3 u(k 1 , m t ), F 12 = u{k 2 , 0)P R u(k 1 , m t ), 
0)P R f(k 3 )^ 3 u(h, m t ), F 14 = u{k 2 , 0)P L f{k 3 )fou(h, m t ), 
0)P R ^*(k 4 )lt 3 u(k 1 , m t ), F w = u(k 2 , 0)P L (/*(k 4 )k' 3 u(k 1 , m t ), 



(2.4) 



The Lorentz invariant t and s are the usual Mandelstam variables, while t 12 = [k\ — k 2 ) 2 and 
s 23 = (k 2 + k 3 ) 2 are the generalized Mandelstam variables. The scalar products are defined 

as 



SP 1 = e*(k 3 )-e*(k 4 ), SP 2 = e*(k 3 )-k 1 , 

SP 3 = e*(k 3 ) ■ k 2 , SP 4 = e*(k 4 ) ■ fci, SP 5 = e*(k 4 ) ■ k 2 . 



(2.5) 



The functions Ri, ...,i?27 are defined in terms of Passarino-Veltman functions |3l| and are 
given explicitly in the appendix. 4 The decay width can be written as 



1 X 

dT{t^cgg) = - — V \M\ 2 d^ 3 {h] k 2 , h, k 4 
2m t z — ' 



spins 



where \Ai\ 2 is straightforward to calculate from the amplitude given in Eq. (|2.3j) . The 
phase space d<&3 can be expressed in terms of energies of the third and fourth particles, 
chosen to be the gluon pair in the rest frame of top quark. In this frame, one can take the 
production plane as the x — z plane and choose the c quark momentum along the z-axis. 
After rearranging the volume elements and carrying out angular and momenta integrals, we 
implement the phase space cuts discussed previously 

d^ 3 (h;k 2 ,k 3 ,h) = — / / dk° 3 I / dk° 2 , (2.7) 



where the limits are 



(k° 2 ) min = Max 



(fc 3 °) min J (*°2) 



Crru.-, — 



2 

(A; 3 °) min = Cm t , 



(^ 3 °) max = ^(1-2C), (2.8) 

with a = m t — k® (recall that for simplicity we have assumed m c = 0). Here C is our cutoff 
parameter, which we initially take as C = 0.001 and then study the effect of increasing its 
value. To calculate the branching ratio, we assume that t — > bW is the dominant decay 
mode of the top quark and use T(t — > bW) = 1.55 GeV. For the numerical analysis, we have 
used the parameters 2^] given in Table ID 



4 The result is expressed in terms of Passarino-Veltman functions [3l| but the reduction to the scalar 
function Aq, Bq,Co, and Dq has not been carried out. This is indeed one advantage of using FormCalc 
which does not require such reduction to save substantial CPU time. 
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TABLE I: The parameters used in the numerical calculation. 



a s (m t ) 


a(m t ) 


sin 0w{mt) 


m c (m t ) 


m b (m t ) 


m t (m t ) 


0.106829 


0.007544 


0.22 


0.63 GeV 


2.85 GeV 


174.3 GeV 



Our result is 

T(t-y egg) 



Br(t -> egg) 



1.02 x 10" 9 , (2.9) 



while for the two-body decay t — > eg we get 

Br(t -> c#) = 5.73 x 10~ 12 . (2.10) 

for the same parameter set. At this point, a comment is needed. For completeness, and to 
check our procedure, we recalculated the two-body decay t — > eg and our numerical value is 
around one order of magnitude smaller than that of the Ref. Q] (see Fig. 2 of The main 
source of such disprepancy lies in the value for bottom quark mass. The pole mass m& = 5 
GeV is used in Ref. jjj, while we have used the running bottom quark mass at m t scale, 
rrif, = 2.85 GeV. Since the branching ratio is proportional to ml (due to GIM suppression), 
our results differ by one order of magnitude ((2.85/5) 4 ~ 0.1). The original SM calculation 
of t — > eg in Ref. jlj was later updated (^J, where the running bottom quark mass was used 
and Br(£ — > eg) = 4.6 x 10~ 12 was obtained. 

There is however one cautionary remark about the t — > egg decay. The branching ratios 
were calculated as function of the cutoff (see below). One should calculate the rate for 
t — > egg by doing a complete higher order calculation. QCD corrections should be at most 
of the order of 10%, which is the order of magnitude of QCD corrections in t — > bW. One 
can view C as a detector cutoff. The cutoff dependence of the branching ratio of t — > egg is 
given in Table ITT1 As seen from Table ITT| the branching ratio is sensitive, but not significantly 
so, to the C parameter. 

In general, contributions from ghost diagrams are quite suppressed with respect to the 
rest of the diagrams and the 1PI diagrams slightly dominate the self energy and vertex type 
diagrams depicted in Figs. [TJ 
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TABLE II: The cutoff dependence of the branching ratio of t — > egg for various C values. 



C (in m t units) 


0.001 


0.003 


0.005 


0.01 


Br(i -» egg) 


1.02 x 10~ 9 


9.04 x 10~ 10 


8.76 x 10~ 10 


7.78 x 10~ 10 



Within the range (0.001 — 0.01) for the cutoff C, the rate for the three-body decay t — > egg 
is more than two orders of magnitude larger than the rate for the two-body decay t eg. 
This is higher order dominance par excellence. However, the branching ratio for the t — > egg 
channel in the SM, is still too small to be observable in any conceivable experiment. Thus 
any experimental sighting of t — > egg will indicate the appearance of NP. The three-body 
decay t — > egg has another important difference with respect to t — ► eg: A promising ratio 
for t — > egg might lead to a sizable cross section for single top production via gg — > tc at a 
future hadron collider, especially since most of the interesting events will be fed by partonic 
sub-processes originating from gluon-gluon collisions. We are currently investigating the 
decay t — > egg and the effect of its crossed partonic sub-proce ss q q — > tc on single top 
production at the LHC, within the minimal supersymmetric SM |3J]. 

III. THE DECAY t -» cqq, q = u 

Using the same procedure and the same parameters as for t — > egg, we have also calculated 
the branching ratio of t — > cqq, q = u. Unlike t — > egg case, this decay arises dominantly from 
the teg* vertex, where g* then decays to — > qq pair. There are of course diagrams mediated 
by electroweak gauge bosons, 7 and Z or the neutral Higgs boson, but their contributions in 
the SM are negligible. The dominant diagrams are given in Fig. El while the box diagrams 
displayed in Fig.@J were found to be at least 3 orders of magnitude smaller than the diagrams 
in Fig. El 

Compared with t —> egg, it is a much simpler calculation, mainly since there are no exter- 
nal gluons and thus no ghosts. The only non-trivial issue is to demonstrate the ultraviolet 
finiteness of the decay. This can be checked either analytically at the amplitude level, or 
numerically for the branching ratio. We followed the second method. 

We find that t — > cqq, q = u remains slightly smaller than t — > eg. Thus we disagree 
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u u u u 

5 6 7 8 

FIG. 4: The sub-dominant one-loop contributions to t — » cqq in the SM in the 't Hooft - Feynman 

gauge. 



with the recent study by Cordero-Cid et al. |20| where t — > cqq, q = u is found to be almost 
an order of magnitude larger than t — > eg. The reason for the discrepancy might be their 
use of different parameters, which unfortunately are not spelled out in their paper. Since we 
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mainly concentrate on the t — > egg decay mode in this study, we prefer to skip the details of 
this calculation and present only our result here. For the parameter set chosen before, we 
found 

Br(t -> cqq, q = u) = 3.96 x 10" 12 . (3.1) 

As seen, unlike Br(t — > egg), this ratio stays slightly smaller but still comparable to Br(t — > 
eg). The decay t cqq has been analyzed and compared with t — > eg in both the SM 
and version II of the Two-Higgs-Doublet model by Deshpande, Margolis and Trottier j^. 5 
They show that for m t = 175 GeV, BR(£ — > cqq) in the SM is slightly (~ 1.2 times) bigger 
than Br(t — > eg) when there is a sum over q = u,c, d, s, b. Our results agree with theirs. We 

5 

get ^2 Br(£ — * cqq) as 1.56 x 10~ n , which becomes bigger than Br(t — > eg). 

9=1 

IV. CONCLUSIONS 

In this study, we have discussed higher order dominance in the rare top decays t — > egg 
within the SM framework, a phenomenon known more than a decade ago in 6-physics. For 
completeness, we also calculated t — > cqq, q = u. Using running quark masses, we have 
found Br(t -> egg) = 1.02 x 10~ 9 with a cutoff C = 0.001 for a top mass m t = 174.3 GeV. 
We considered the sensitivity of the ratio with respect to the cutoff parameter C and found 
that even for C = 0.01 it is still more than two orders of magnitude larger than the two-body 
decay Br(£ — > eg) which we calculated as 5.73 x 10~ 12 with the same set of parameters. By 
comparison, we found Br(£ — > cqq, q — u) to be smaller, 3.96 x 10~ 12 and comparable to 
Br(£ — > eg). However, when we sum over qq pairs for all five quarks, it becomes slightly 
larger than Br(£ — > eg). 

If higher order dominance is still valid for a viable NP model in the sense that t — > eg is 
much smaller than t — > egg yet larger than its value in the SM, we may have a glimpse of 
NP at work either in the decay t — > egg or in production at the LHC through the partonic 
sub-process gg — > tc. 

5 Some of this decay mode (q = d, b) has also been considered by Eilam et. al. in the SM [3f]|. 
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APPENDIX: THE ONE-LOOP R-FUNCTIONS 

The explicit form of the functions defined in Eq. (|2.3|) are given below. As in the paper, 
expressions are given for m c = 0, while we use a non-zero c-quark mass in our numerical 
calculations. We define T d = (T a T^)ij, T e = (T\,T a )ij as products of color matrix elements. 
Here a and b are gluon indexes running from 1 to 8 with color indexes i,j = 1, 2, 3 

Ri = m 2 t - s 23 , 

R 2 = 2m 2 t m 2 w s 23 (-tT d + m 2 t T e - s 23 T e ) + ((ml - m^)(m& + 2m 2 v )(m 2 - s 23 )s 23 T e 
-m 2 t (m\ - 2m% / + m 2 b (m 2 w - 2s 23 ))tT d )B^ 1] - (ml + 2(m 2 - m 2 w )m 2 w 
+m 2 b (m 2 v - m 2 ))(m 2 - s 23 )s 23 T e Bf ] + 2m 2 b m 2 s 23 (s 23 - m 2 t )T e Bf> 
+m 2 t (ml + m\(m 2 w - s 23 ) + 2m 2 v (s 23 - m^))tT d B^\ 

R 3 = (m\ - m 2 w )(m 2 b + 2m 2 w )(tT d - (t + 2t 12 )T e )(E< 1 1) - B {2) ) + m 2 t (m 2 b (tT d (B {2) 
-2B { 1] ) + T e (2tB { 1] - tBf - 2t n B® + 4t n B®)) 
-2m 2 w (tT d - (t + 2t 12 )T e )(B {2) - 1)), 

Ri = m 2 t — t, 

R 5 = 2m\m 2 w ((m 2 t - t)t(T d - T e ) - (m 2 t - 2t)t l2 T e ) + (m 4 b + m 2 b (m 2 w - 2m 2 ) 
-2m 4 w )t((m 2 t - t)(T d - T e ) + t 12 T e )B^ - (m 4 b + 2(m 2 t - m 2 w )m 2 w 
+m 2 b (m 2 w - m 2 t ))(m 2 t - t)(tT d - (t + t 12 )T e )B {2) + m 2 t (-m\ - 2m 2 w (t - m 2 w ) 
-m 2 b (-2m 2 + m 2 w + t))t l2 T e Bf\ 

R 6 = 2m 2 w (m 2 (m 2 t - 2s 23 )t 12 T d + m 2 (m 2 t - s 23 )s 23 (T d - T e )) + (ml + m 2 b (m 2 w - 2m 2 ) 
-2m 4 v )s 2 3((m 2 - s 23 )(T d - T e ) - t 12 T d )B^> - (ml + 2(m 2 t - m^)m^ 
+ml(m 2 v - m 2 ))(m 2 - s 23 )(t l2 T d + s 23 (T d - T e ))B^ 2) + m\(ml 
+2m 2 w (s 23 - m 2 w ) + m 2 b (-2m 2 t + m 2 w + s 23 ))t l2 T d Bf\ 

R 7 = 2m 2 m% / (m\(tT d + s 23 T e ) + 2s 23 t(tT d + s 23 T e ) - m 2 t (T e s 2 23 + 2t(T d + T e )s 23 
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+t 2 T d )) - (m 4 b + m 2 {m 2 w - 2m 2 t ) - 2m 4 w )s 23 t(-tT d - s 23 T e + m 2 {T d + T e ))B^ 
-{ml + 2(m 2 - m\)m\ + m^ h {m^ w - m 2 t )){m 2 t - s 23 )(m t 2 - t)(tT d + s 23 T e )B^ ] 
+m 2 (m 2 - s 23 )s 23 {m 4 b + 2m^(t - m 2 w ) + m 2 b (~2m 2 t + m 2 w + t))T e Bf> + m^(m£ 
+2m^(s 23 - m^) + m 2 b (-2m 2 + m 2 w + s 23 ))(w^ - t)tT d B^\ 

2m 2 m 2 w (s 23 (s 23 - m 2 )(SP 2 - SP 3 )t(t - m 2 )(T d - T e ) - t 12 {m 4 t {SP 3 tT d + s 23 SP 2 T e ) 

+2s 23 t(SP 3 tT d + s 23 SP 2 T e ) - m 2 (SP 3 t(2s 23 + t)T d + s 23 SP 2 (s 23 + 2t)T e ))) 

+ K + m 2 b (m 2 w - 2m 2 ) - 2m 4 w )s 23 t((m 2 - s 23 )(SP 2 - SP 3 )(m 2 - t)(T d - T e ) 

+t 12 (m 2 SP 3 T d - SP 3 tT d + m 2 SP 2 T e - s 23 SP 2 T e ))Bl 1] - (ml + 2(m 2 - m 2 w )m 2 w 

+m 2 b (m 2 w - m 2 ))(m 2 t - s 23 )(m 2 - t)(s 23 (SP 2 - SP 3 )t(T d - T e ) - t 12 (SP 3 tT d 

+s 23 SP 2 T e ))B { 2) - m 2 (m 2 - s 23 )s 23 SP 2 (m 4 b + 2m 2 v (t - m 2 w ) 

+m 2 b {-2m 2 t + m 2 w + t))t 12 T e B { 3) - m 2 t im\ + 2m 2 v (s 23 - m 2 w ) 

+m 2 b (-2m 2 + m 2 w + s 23 ))SP 3 (m 2 - t)tt 12 T d B^ , 

s 23 SP,T e {2m 2 w {C? + tf«) - rn^cf + Cf - C§)) + SP,(m 2 - t) 

xT d (2m 2 w (C$ + C%) + m 2 b (C^ + C$ + 2^ 4) + eg)), 

SP 3 tT d {m 2 b {Cf ] + cf - cg>) - 2m 2 w (C{ 3) + C@>)) - (m 2 - s 23 ) 

xSP 2 T e (2m 2 w (cW + Cif) + m^cf + + 2Cf + Cf)), 

m 2 w {m 2 - t)t(T d - T e ) - (m 2 + 2m 2 w ){m 2 - t)tB^{T d - T e ) 

+ (m 2 - t)t{{{m 2 b -m 2 - m 2 w ){m 2 b + 2m 2 w ) + 2m 2 w t l2 )C^ + 2{m\ + 2m 2 w )Cf 

+2m 2 w (t 12 - m 2 t )Cf ] - {m 2 b m\ + Am 2 w m 2 - 2m 2 w t 12 )C^)(T d - T e ) 

-m 2 w (m 2 - 2t)t 12 T e + (m 2 b + 2m 2 w )(m 2 - 2t)t l2 T e Bf ] + tt 12 T e (-m 2 b m 2 C^ 

+2(m 2 + 2m 2 w )C ( £ ) + 2m 2 w (m 2 - s 23 - t 12 )C?) + (m 2 - t)t l2 T e (m 2 m 2 C^ 

-2{m\ + 2m 2 w )Cf ) + (m 2 b m 2 - 2m 2 w t)C { 2 2) ), 

m 2 w (-tT d + m 2 T e - s 23 T e ) - (m 2 b + 2m 2 v )(m 2 - s 23 )T e B^ } + (m 2 b + 2m 2 v ) 
xtT d B® + m 2 b [m\ - s 23 )tT e c£ ] + tT d {m 2 b s 23 C^ - 2(m 2 b + 2m^)Cj 3) 
-2m 2 w s 23 Ci 3) ) - (m 2 - s 23 )T e (m 2 b ((m 2 - t)(Cf } + Cf ) - 2Cf ) 
-4m 2 w Cf J) - (m 2 - s 2 - 3 )tT d {m 2 b {C^ + ^ 4) ) - 2m 2 w C^), 
m 2 v s 2 3(s 23 - m 2 t )(SP 2 - SP 3 )t(t - m 2 t )(T d - T e ) - (m 2 b + 2m 2 w )(m 2 t - s 23 ) 
xs 23 (SP 2 - SP 3 )(m 2 - t)tBf\T d - T e ) + (m 2 - s 23 )s 23 (SP 2 - SP 3 )(m 2 - t)t 
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x(((m 2 - m 2 - m 2 w ){m 2 + 2m 2 w ) + 2m 2 w t 12 )Cf ] + 2(m 2 + 2m 2 w )Cf 

+2m 2 w (t 12 - m 2 )C[ 5) - {m 2 b m\ + Am 2 w m 2 - 2m 2 w t 12 )C^){T d - T e ) 

-m 2 w t 12 (mt(SP 3 tT d + s 23 SP 2 T e ) + 2s 23 t(SP 3 tT d + s 23 SP 2 T e ) - m 2 t (SP 3 t(2s 23 + t)T d 

+s 23 SP 2 {s 23 + 2t)T e )) + (m 2 + 2m 2 w )K - s 23 ) S23 SP 2 ( mt 2 - 2t)t 12 T e B^ ] 

+ {m 2 + 2m 2 w )(m 2 - 2s 23 )SP 3 {m 2 - t)tt 12 T d B® - {m\ - s 23 )s 23 SP 2 tt 12 T e 

x {m 2 b m 2 t C^ - 2{m\ + 2m 2 w )C^ + 2m 2 w (-m 2 + s 23 + t 12 )C[ 1] ) 

+s 23 SP 3 {m 2 - t)tt 12 T d {2{m 2 b + 2m 2 w )Cf ] + m 2 b m\cf - (m 2 + 2m 2 w )s 23 C[f) 

+ (m 2 - s 23 )SP 3 (m 2 - t)tt 12 T d (m 2 b m 2 C^ - 2{m\ + 2m 2 w )Cf + (m 2 m 2 

+2m 2 w (-m 2 + t + t l2 ))C { 2 ] ) - (m 2 - s 23 )s 23 SP 2 (m 2 - t)t 12 T e 

x (2(m 2 + 2m 2 w )C {2) + (m 2 + 2m 2 w ){m 2 - t)C[f + m 2 b m 2 t C { 2 2) 

+m 2 t {m 2 b + 2m 2 w )Cf 2 ) ), 

m 2 w {m\{tT d + s 23 T e ) + 2s 23 t(tT d + s 23 T e ) - m 2 t {T e s 2 23 

+2t{T d + T e )s 23 + t 2 T d )) - {m 2 h + 2m 2 w ){m 2 - s 23 )s 23 {m 2 - 2t)T e B^ 

-{m 2 b + 2m 2 w ){m 2 - 2s 23 ){m 2 - t)tT d B^ ] + {m\ - s 23 )s 23 tT e {m 2 b m 2 C^ 

-2{m 2 b + 2m 2 w )C^ + 2m 2 w {-m 2 + s 23 + t 12 )C[ 1] ) + s 23 (m 2 - t)tT,(m 2 m 2 ^ 3) 

-2(m 2 + 2m 2 y )C,S 3) - 2m 2 w s 23 C[ 3) ) - {m\ - s 23 )s 23 {m\ - t)T e {m 2 b m 2 t C {2) 

-2(m 2 + 2m 2 w )C {2) + {m 2 b m 2 - 2m 2 y t)C? } ) - {m 2 - s 23 ){m 2 - t)tT d 

x {mlmlcf - 2(m 2 + 2m 2 w )Cf + (m 2 m 2 + 2m 2 y (-m 2 + t + t 12 ))C^), 

m 2 w {SP A - SP 5 ){{m 2 - 2s 23 )t 12 T d + {m 2 - s 23 )s 23 (T d - T e )) - (m 2 + 2m 2 ,) 

x(m 2 - s 23 )s 23 (SP 4 - SP 5 )(T d - T e )B^ - (m 2 + 2m 2 v )(m 2 - 2s 23 )(SP 4 - SP 5 ) 

xt 12 T d B^ + s 23 {SP 4 - SP,)t l2 T d {m 2 m 2 C^ - 2{m 2 + 2m 2 v )Cf 

-2m 2 w s 23 C[ 3) ) - {m 2 - S2 3)s 2 35P 5 ti2T e (m 2 ci 2 ) + 2m 2 y (C , { 1) + Cg } )) 

+ K - s 23 )t 12 r,(-m 2 m 2 ( 1 SP 4 - ^P 5 )Cj 4) + (m 2 + 2m 2 y )(2(,SP 4 - 5P 5 )Ci 4) 

-s 23 ^P4Ci 4) ) - {m 2 b m 2 {SP 4 - SP 5 ) + 2m 2 w (m 2 (SP 5 - SP 4 ) - SP 5 {t + t 12 ) 

+SP 4 (s 23 + t + t l2 )))C { 2 A) ) + {m 2 - s 23 )s 23 (SP 4 - SP 5 )(T d - T e ) 

x {{{m 2 b - m 2 - m 2 w ){m 2 b + 2m 2 w ) + 2m 2 y t 12 )C'J 5) + 2(m 2 + 2m 2 w )C^ 

+2m 2 w {t 12 - m 2 )C[ 5) - {m 2 b m\ + 4m 2 y m 2 - 2m 2 w t 12 )Cf ), 

m^K 2 - s 23 )s 23 SP l {m 2 - t){tT d - {t + 2t 12 )T e ) + 2(m 2 + 2m 2 w ) 
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x(m 2 - s 23 )s 23 SPi(m 2 - t)t 12 T e B { 3) - (ml + 2m 2 v )(m 2 - s 23 ) 
xs 23 SP 1 (m 2 t - t)t(T d - P e )pf - 2(m 2 t - s 23 )s 23 tt 12 T e (m 2 b (m 2 SP 1 Ci j 1) 
-SptCf + 2SP 2 SP 5 (C { 1) + Cf - C$)) - Am 2 w SP 2 SP,(C? 
+C$)) - As 23 SP 3 (SP 4 - SP 5 )(m 2 - t)tt 12 T d {m 2 (C^ + cf - Cf) 
-2m 2 w (C[ 3) + Cf)) + 2(m 2 - s 23 )s 23 SP x {m 2 t - t)t 12 T e (m 2 ((m 2 - t)(C< 2 > 
+cf ) - 2Cf ) - 4m^Cf ) + 4(m t 2 - s 23 )SP 3 SP 4 (m 2 - t)tt 12 T d (2m 2 w (C$ 
+Cff) + ™ 2 (Cf + Ci ( ? + 2C 2 (4) + C*f)) + (m 2 - S23 )s 23 (m 2 - t)t(T d - T e ) 
x{m i b SP 1 C§ ) + 2m 2 w (-m 2 w SP 1 C i i ) - 2s 23 SP l C { * ) - ASP 2 SP b C§ ] + 2SP 1 Cf ) 
-2s 23 SP 1 Cf ) - 4 1 SP 21 SP 5 C , { 5) + rrfiSPxC^ - 2s 23 SP 1 Cf 2 ) - ASP 2 SP^Cf 2 ] 
-SPit 12 C$ - 2m 2 t SP 1 C { 2 ) - 2s 23 SP 1 C< 2 5) - 8SP 2 SP 5 Cf } + 2SP 1 tC { 2 ) + SPA^f 
-(ASP 2 SP 5 + SP 1 (m 2 -2t- t 12 ))C$ + ASP 3 SP 4 (C® + c[ 5) + c[f + 2Cf } 
+Cf )) + m 2 b {A{SP 3 SP, - SP 2 SP,){-Cf ] + Cf + Cf ) - m 2 SP 1 (Cj ) 5) + cf 
_tfg> + Cf + cf) + SP^Cf + 2Cf + (2s 23 + t 12 )(Cf } - cg>) 
+(2t + t 12 )cg ) ))), 

2m 2 w SP 1 T d {Cf ) - cf ) + mlSPMC? > - Cf) + 2T d (m 2 (5P 1 (p( 1) + Pf 
+m 2 w (D? - Pf )) + SP 3 (SP 5 (Pf 2 + Pf + Pg 2 ) + SP 4 (Pf 2 + Dg - Pf 
+Pf + Pf 2 + Pf 3 )) + ^^(^(DS + Pf ) + SP 4 (Pf 3 + Pf + pf 3 ))) 
+mym^P 1 (P[ 1) - Pf ) + 2(SP 3 (£P 5 (Pf 2 + £>f + D$ 2 ) + SP 4 (2Pf 
+#£a + D<& + Pf + 2Pf + 3Pf + pf 2 + pf 3 + 2Pf )) + SP 2 (SP 5 (pf 3 
-Pf ) + SP 4 (Pf 3 + Pf + P£ 3 ))) + SP^-W™ + 2Pf - S23 (pf + 2Pf ) 
+ (-2m 2 + t + i 12 )pf + m 2 Pf ))) + TMSP^D® - Pf ) - m^m^A 
x(pf - pf + pf - pf ) + S'P 1 (-4PS 2) + 2Pf + (s 23 + t 12 )pf + m 2 w {Df ] 
-Pf ) - i(Pf + 2Pf )) + 2(SP 3 (SP 4 (Pf + Pf 3 ) + SP 5 (P? 13 + Pf 
+P{|)) + SP a (SPt(D® + Pg + Pg 3 ) + ^(-M? + + + Df,\ + Pg 
+Pg 3 )))) - 2m 2 y (m 2 5P 1 (Pf + Pf } + Pf + Pf ) - SP(t 12 (Pi 2) + Pf } 



f + Pf ) - 2(Pf + Pf )) + 2(SP 3 (SP 4 (Pf 23 - Pf ) + SP 5 (Pf 3 + Pf 
f 3 )) + 5P 2 (5P 4 (Pf 3 + Pf + Pf 3 ) + 5P 5 (Pf + Pf 3 + 2Pf 

(2) , on (2) , n(2) , on(2) , o r»(2) , r»P) 



-pr 33 + 2P^ + p 233 + 2Pf + spf + pf 3 ))))), 
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m 2 (SP 3 (T e (D$ - D?) + T d (-D^ + D$ + flg + D% + £>£>)) 
+SP 2 (T e (D$ + Dfg - D®) + T d (D® + -D33 ))) + 2m^(5P 3 (T d (4 1) + D? J 

- D« - D% - Dl? + D«) + UDf + ) - 2>g> + + D< 2 >)) 
-SP 2 (T e (D? + D® + £>£>) + T d (D® + L>« + £©))), 

m^(T d (D« - D«) + T B (D< 2) - + m^cf - C? - 6D^ + 3m^(I> W 

- ^(Dg + 2D$) + (t + t 12 )D$ + m?(U w - 2D$ + 
+T e (Cf - cf - 6L>f + (s 23 + toa)^ + m 2 ^ - D® 
+D® - Dg>) + 3m 2 w (D^ - D® ) - f(Z>g + 2Z>g>))) + 2m^(T d 
x(-Cf + Cf + 6Djp + - + s 23 (dS + 2I>£>) 

+ ml[D^ + 0« + + 21$) - i u (!)« + D? + + D$ 
+D®)) + U-CP + cf + 6ltf> - ajB i?g) + m 2 w (D^ - D®) - MD® + D$ 2) 
+-D13 + + D?) + m 2 ^ + D? + D% + D% + D$>) + t(D$ + 2Dg>))), 
2m 2 w {SP,(T d {D^ + ^ - D$ + D« + D s «) + r e (D « + - D® - D<$ 

- - £>g>)) - SP^UDf + Z>g> + £>g>) + T d (L»« + D 8 W + £>«))) 
+m 2 (5P 4 (T e ( J Dg ) + Z>g>) + T d (-D« + -D23 + ^3)) + ^WDff - 
+T e (-Df) + Dg) + Pf 3 ) +^ 2 3 ) +^ 2 3 ) ))), 

T d {2m 2 w {C^ + Cf ) - m 2 Cf) + T e ((m 2 - 2m^)(C{ 7) + cf) - m 2 Cf ) 
+T d (m^« + 2m^(-2D« - s 23 (£>« + L>« - D« - 21$) + m 2 t D^ 
+m 2 w {D ( i ) + £>«) - (-2m 2 + t + f 12 )D«) + m 2 (2Z$ + s^D? - + 
-21$) - m 2 w (2D^ + 3Z$) + (-2m 2 + t + t 12 )r$)) + T e (mi(D^ - £>f 

) + m 2 (4Z$ + ( S23 + ^(fl? - D® - D$) + m>(-D® - D? + 2?g 
^f) - Pf + + 0g?) + n4(3(U a W + D$>) - D^) + t(I?g) + 3D$ - D$» 
+2DS)) - 2m^(4Lf - t 12 (L>S? + £>g>) + m 2 (L>i? + D§ + £>g> + D§) 
+m 2 w (D® + 2?3 W ) + «23(^ 2) + D? - D® - D® + + ) + t(Df + 2?g> 
+3Dg ) + D a ca) + 2Dg ) ))) ) 

2mt v SP 1 T d (Di 1) + D®) + m 2 b (-SP 1 (2T d (Di 1) - Z>«) + T e (Cf } - - Cf 
+2(-2D^ + + Df) - ( S23 + - - D$) + m?p « + p! 2) 

-D® - D% + D? - D$ - D$) + m 2 (- D W + D W + D f) _ t{D $ 
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+3Dg> - D® + 2D<§))) - 2(SP 3 (SP 5 (T e (-D ( $ + + D W + D M + Dfi) 
-T d (L>« - < + 2£>« - < + + + SP,{T e {D?l + L>? 23 ) 

+T d (- J D« + + + Df 3 » - + !>£> + 2I>£> + + D<?))) 
+SP 2 (SP,(T e (D^ + + 2£# 3 + £>g> + Dg,) + UDg\ + + £>§!,)) 
+SP 5 (T e (-DP + £>§> + 2£>g 3 + Dj» + Z)g> + 2?g> + 2^1 + £>g> + D«) 
-T d (I$ - + D% + £>«))))) + m^SP^UC^ + Cf - 2D^ - 2L>« 
+s 23 (dS + + 2D%) - ml(D^ + D™) + (2m,? - t - + 
-T e (4(itf> + + D?) + ml(D^ - D® - D?) + 2s 23 (dP + D« + 
+D®))) - 4{SP a (SP B (T e (D® + D$ + + P{? 3 + 2Z>g> + £>? 23 + 2D (2) 
+Z>£) + T d (P« + D$ + Z?g> + 2 J D? 2 ) + D® + 2^ + I#> + Z>£> + 
+SP A (T e (D® + D§> ) + T d (L>« + D« + + D« + J3« + + 
+2£>« + + 2DP + £>«))) + SP 2 {SP A (T e {D ( £ + D® + 2Z>g> + D® + 
+T d (L»S + D§> + I>gj> )) + SP 5 (T d (D$ + + D« + D<p + £>g) + T e (D® 
+D{ 2) + 2£>{ 2) + Dg + 2£>g>, + 2P© + Dfl + 3D< 2) + 2£>g } + L>g 3 
+5Z>g> + 2Z>g> + 3(£>< 2) + Dg) + I>g> ))))), 

2m' w (SP 2 (T d DP - T e (D { p + D® + Pf)) - SP 3 (T d (D^ + D< x) + D^) 
-T e D { p)) + m2(-5P 2 T e Ci 7) - SP 3 {T d {C^ - m 2 ^ - 4D« - 2D® - 2D? 
-s 23 (dS + D$) + (t + tu)(D<8 + D«) + ml{DP + !>« - dJ? + D« - Z?g> 

+dP)) - u&p + Cf + cf) + 2Df) + K - - t 2 )^J 2) - t(^f 2 } + D?h)) 

-SP 2 (-T e (m 2 b D { p + 2(2Pf + D® + ) - (s 23 + f 12 )Df 2 } + m 2 t (-D^ 
+Dfi + D$)) + T d (C® - 2D? + s 23 (D<$ + D«) - + m^p + 

+T e ((, 23 + t^ug + m 2 (M 2) - - D« + L»f ) - t(D% + 
K (5P 3 (T d (2Cf + ml(DP + 2(P« + DP)) + 2(2(PS 1) + D«) - s 23 (D$ 
+ DP + 2P«) + + I?W + DP + P«) + irf + 2P« + D{? + 2DP 

+ DP + 2P« + Dg> + D 8 «) - m 2 (P« + 2DP + Dj? + 3D« + + 3P« 
+P« + 3P« + 2D?>))) + T e (2Cf - m 2 (2Pf + D?) + 2(-2Pf + 2D? } 
+s 23 D$ + (t + t l2 )(D% + Dg?) + tPf + m 2 (P^ + DP - DP - DP + 2DP 
+DP)))) + SP 2 (-2T e (cS ] + C[ 7) + Cf } ) + T d (2C { P - m 2 b (2D^ + £>«) 
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+2(-2D? + 2D? + s 23 D { V + t 12 D$ + s 23 D? + (s 23 + t + t 12 )D? - m 2 (D^ 
+2D? + D?))) + T e (m 2 (D? + 2{D? + D?)) + 2(2(L>f + L>f ) + s 23 (D^ 
+D? + I?g> + Df) - m t 2 (Df + D? + D? + i?g> + ™? + 2I?g> + 3D? 
+2D?) + t 12 (D< 2 > + D? + + Dg> + 2D? + D« + 2(D? + D< 2 >)) 
-HDtg+D® + 21$))))), 

-m 4 b (SP,(T d D? + T e Z>f ) + SP 5 (T d (-D? + I)' 1 ) + D«) + T e D?)) 
+2m 2 w (SP 5 (T e (C? - C? - cf - 2(3D< 2) + L»J 2) + D?) + m 2 w/ (2^ 2) + D? } ) 
-t 12 (Df ) + D<?) + S2 3^f 3 } + -?(^ 2) + 2^i 2) + 2>g> + D? + 3D< 2) + ^ 2) ) 
+*pg> + 2(2>g> + Pf))) - Ztfcf + 4D« + 2D? + + S23 (D« + i#> 

+2D«) - t(Z>g> + »?) - t 12 {D? + Uf) + Z>« + ^ + D?) + m?(J#> 
+-Di^ + D? + L>« + 2D? + D^))) + SP 4 (T,(Cf } + C? + C? 5 - 2D? 
-2D? + Sa3 (D« + D? + 2(D« + D? + !?«)) + m 2 ,^ + D? + D?) 
+m 2 (D? + D? + D« + ^ + ^ + Ug } ) " tn(D? + 2>?> + D« 
+D« + 2D? > + D®)) - T e (C? + m 2 w D? - 2D? + 2D? - s 23 D$ + m 2 D? 
+tD?))) + m 2 b (SP,(T e (C? - 2D? - (s 23 + t 12 )D? + m 2 w (2D? 
+ D?) + tD? + m?(Dg> + D? + 2I?g + Dg>)) + T d (C? - AD? - 2D? 
-2D? - s 23 {D? + D? + flg> + D«) + m 2 (D? + D?) - m 2 w {D? + 2(D? 
+D?)))) + SP 5 (T e (C? + 2D? - 2D? - 2D? + m 2 w (-AD? - 2D? + D?) 
+t(D? - D?) + m 2 (D? + I?g> + 2D? + ^ + Z)g>)) - r d (2Cf + Cf - 6D? 
+2D? + m 2 w (D? - D? - 3D?) - 2s 23 {D? + D?) + (t + t 12 )(^ 1 3 ) 

+^23) + m t( D ? - 2 ( D S + ^S) + D ?))))i 

Td — T e , 

m 2 b {s 23 {C? - C?) + (m 2 - t)C?) + 2m 2 w (s 23 (C? + C? + C? + C?) 
+ (m 2 -t)(C? + C?)), 

{ml + merrily - 2m 2 ) - 2m^)B? - (ml + 2{m\ - m 2 w )m 2 w + m 2 h {m 2 w - m 2 t ))B? 
+2m 2 {m 2 b + 2m 2 w )B? + 2m 2 (-{(m 2 b - m 2 - m 2 w )(m 2 b + 2m 2 w ) + 2m 2 w t 12 )C? 
-2(m 2 b + 2m 2 w )C? + 2m 2 w (m 2 - t 12 )C? + (m 2 b m 2 + Am 2 w m 2 - 2m 2 w t 12 )C?). 
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where 



A 



B {m 2 b , m 2 w ), = B (m 2 , m 2 b , m 2 w ), B^> = B (t, m 2 b , m 2 w 



>(3) 



B (t 12 ,mlmi), B^ = B (s 23 ,m 2 b ,m 



,2 ™2 



^(0, t, 0, m 2 b , m 2 w , m 2 b ), Cf ] = C P {Q, t, m 2 , m 2 b , m 2 b , m 2 w ), 
Cp(0, s 23 , 0, rri 2 , m 2 w , m 2 b ), Cf = Cp(0, s 23 , m? > m &> m &> m u0 ; 
(7/3(0, iu, m 2 , m^, m?, m 2 ), C^ 6) = Cp(0, 0, ti 2 , m^, m 2 b , m 2 b ), 
C p {t l2 , 0, 0, m 2 w , m 2 b , m 2 b ), (3 = 0, 1, 12, 2, 22, 
L> A (0, 0, 0, m 2 , s 23 , *i2, m^, ml, m 2 b , ml), 
D x (0, t 12 , 0, t, m\, 0, m^, m 2 b , m 2 b , m 2 b ), 

0, 1, 11, 12, 13, 2, 22, 23, 3, 33, 003, 112, 113, 122, 123, 133, 222, 223, 233, 333. 



The coefficient functions C^- and D^ij^jh are symmetric functions and can further be de- 
composed into the scalar functions A ,B ,Co, and D . See j3| for details. 
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